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Abstract 

Exact dynamical equations for a generic dust matter source field in a cosmological context are 
formulated with respect to a non-comoving Newtonian-like timelike reference congruence and investi- 
gated for internal consistency. On the basis of a lapse function N (the relativistic acceleration scalar 
potential) which evolves along the reference congruence according to TV = a O TV (a — const) , we find 
that consistency of the quasi-Newtonian dynamical equations is not attained at the first derivative 
level. We then proceed to show that a self-consistent set can be obtained by linearising the dynamical 
equations about a (non-comoving) FLRW background. In this case, on properly accounting for the 
first-order momentum density relating to the non-relativistic peculiar motion of the matter, addi- 
tional source terms arise in the evolution and constraint equations describing small-amplitude energy 
density fluctuations that do not appear in similar gravitational instability scenarios in the standard 
literature. 
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1 Introduction 

The aim of this paper is to investigate the Newtonian limit of the Theory of General Relativity in a 
cosmological context, when conditions relating to gravitational physics are significantly different from 
the more traditional quasi-stationary and asymptotically flat situations where the Newtonian and post- 
Newtonian limits can usually be derived and are subject to precise experimental tests. The importance of 
this is that many astrophysical calculations on the formation of large-scale structure in the Universe rely 
on such a limit, because most of these calculations are done in a Newtonian way. A number of derivations 
of such a limit are given in the literature, see, e.g., Refs. j| and (^7); indeed, there is a viewpoint 

that cosmology is essentially a Newtonian affair, with the relativistic theory only needed for examination 
of some observational relations. However: 

(A) The true classical theory of gravitation is General Relativity; Newtonian theory is only a good 
theory of gravitation when it is a good approximation to the results obtained from General Relativity; 

(B) We have to extend standard Newtonian theory (which strictly can only deal with quasi-stationary 
isolated systems in an asymptotically flat spacetime) in some way or another to deal with a non-stationary 
cosmological context. This extension needs to be clearly spelt out. 

Significant questions remain in regard to both points. These include: 

(1) Full General Relativity theory involves ten gravitational potentials (combined in a tensorial vari- 
able), subject to the ten Einstein field equations ('EFE'), but Newtonian theory involves only one (a 
scalar variable), subject to the one Poisson field equation ('PFE'); how does it arise that the other nine 
potentials and equations can be ignored in the Newtonian limit? Given that nine equations of the full 
theory are not satisfied even in some limiting sense, how do we know when Newtonian cosmological 
solutions correspond to consistent relativistic solutions of the full set of equations? Part of the answer 
is that the coordinate freedom of General Relativity accounts for four of these potentials; however, this 
still leaves another five to account for, and we have examples of Newtonian solutions with no General 
Relativity analogues |2(J . Consequently we need to be concerned as to how well standard Newtonian 
theory represents the results of General Relativity in the cosmological context we have in mind, when we 
take a 'Newtonian limit'. 

(2) The issue of boundary conditions for Newtonian theory in the cosmological context is problematic 
even in the context of exactly spatially homogeneous (and spatially isotropic) cosmological models J26| ; 
no fully adequate theory exists in the more realistic almost-homogeneous case. Numerical simulations, for 
example, usually rely on periodic boundary conditions, which correspond to the real Universe only if we 
live in a 'small universe' |]l6| in which there is a long-wavelength cut-off in the spectrum of inhomogeneities 
of its large-scale structure. Analytic solutions usually rely on asymptotically flat conditions which are 
manifestly not true in a realistic, almost-Friedmann-Lemaitre-Robertson- Walker ('FLRW') situation 
(they implicitly or explicitly assume that inhomogeneities die off sufficiently far away from the region of 
interest). 

(3) How do we attain a unique propagation equation for the gravitational scalar potential in a New- 
tonian cosmology, when Newtonian theory proper has no such equation? In standard Newtonian theory, 
which aims to describe only quasi-stationary settings, this results from the assumed asymptotically flat 
condition — which does not hold in the cosmological context. 

(4) How do we satisfactorily handle the gauge- dependence that underlies most derivations of a New- 
tonian limit? Equivalently, most derivations of the Newtonian limit are highly coordinate dependent, 
despite the strongly proclaimed doctrine of covariance of General Relativity; do we need to revise our 
view on this?f] 

One approach to these issues is the 'frame theory' suggested by Ehlers [|[ ||, which has the advantage 
of giving precise theorems on some Newtonian limits, but has the disadvantage of being more complex 
than General Relativity. Here we seek a more direct approach to the Newtonian limit in a cosmological 
situation, that is more readily related to the approaches taken in the astrophysical context. We do so by 
introducing a class of exact relativistic cosmologies that we can justifiably call quasi- Newtonian because 
of the properties we discuss below, and which include as special cases linearised models used in many 
astrophysical studies, for example in gravitational lensing theoryj^j In studying these spacetimes, we will 

1 Ideally one would impose conditions at a finite distance rather than infinity jl3| ; however, this 'finite infinity' approach 
has not yet been developed fully. 

2 See, e.g., the discussion by Ellis and Matravers 1191 . 

3 They are related to the spacetimes given this name by Bardeen [p| when we linearise our models about a non-comoving 
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come across each of the above issues, in a context where we can pose — but not often solve — precise 
questions concerning exact quasi-Newtonian solutions of the EFE; consequently we will present some 
unanswered 'Problems ' in what follows — well-posed questions that should be answerable, but are of 
such complexity that we have been able to make only little progress in their solution. We then examine 
the effect of linearising about a (non-comoving) FLRW model, when these questions become tractable 
and interesting issues occur in relation to the approaches usually taken. 

In order to limit the complexity of the problem, we will only consider cosmological models in which the 
source of spacetime curvature is given by a pressure-free matter field ('dust'); that is, we will examine here 
only the pure gravitational problem. A further feature of our discussion will be the absence of gravitational 
radiation; on the one hand, this is a property of gravitational fields unknown to the Newtonian theory, 
and on the other, gravitational radiation is commonly thought of as not being of immediate dynamical 
relevance in structure formation processes. We proceed as follows: in section || we develop our concept 
of a quasi-Newtonian dust cosmology in 1 + 3 covariant terms and touch upon issues of the uniqueness 
of the non-comoving time slicing chosen. In section [| we present the relevant 1 + 3 covariant dynamical 
equations. These are derived from the Ricci and Bianchi identities on using the EFE to algebraically relate 
the Ricci curvature to the dust matter source field. We then fix the time-reparameterisation freedom 
of the EFE by proposing to evolve the lapse function along the timelike reference congruence in a way 
which includes various choices made in the literature. In section^, with the given choice of lapse function 
evolution, we pursue a consistency analysis of the dynamical equations to the first derivative level; no 
consistency is attained at this stage. A number of dynamically specialised subcases are then highlighted 
in section ^[ In section ^| we investigate the effect of consistently linearising our dynamical equations 
about a non-comoving FLRW background. Compared to results in the standard literature, our approach, 
which fully accounts for the first-order momentum density of non-relativistically moving matter, leads to 
modifications in the set of evolution and constraint equations. Our conclusions are contained in section ^. 
Useful computational details such as commutation relations, the set of 1+3 covariant dynamical equations 
cast into an associated 1 + 3 orthonormal frame form, and the linearised EFE for a scalar-perturbed, 
spatially flat FLRW model have been gathered in an appendix. Throughout our work wc will employ 
geometrised units characterised by c = 1 = 8nG/c 2 . Some of our notation is inspired by Ref. |3C[ |. 

Our discussion concentrates on two main issues: (i) the internal consistency of the set of exact equa- 
tions defining relativistic quasi-Newtonian dust cosmologies; in looking at this, we need to ask whether 
or not the time slicing condition we impose leads to a uniquely defined foliation T: {t — const}, and have 
to face the problem of time-reparameterisation freedom of the dynamical equations of General Relativity 
(the EFE); and (ii) how, if at all, do the exact equations reduce under FLRW- linearisation to the New- 
tonian form encountered in many works on gravitational instability of matter in non-relativistic motion, 
i.e., what happens to the constraints in this limit, and when is the set of equations so obtained both 
consistent and linearisation-stable? Of most direct relevance to an astrophysically oriented reader is the 
discussion given in section ^, where, in relating our approach to some of those given in the astrophysical 
literature, we show that extra terms not always included seem to be needed if the Newtonian-like equa- 
tions used are to be consistently derived from General Relativity. The reader may directly address this 
material after becoming acquainted with the definitions and relations provided in sections and 0. 



2 Quasi-Newtonian cosmologies 
2.1 Definition 

We define relativistic quasi- Newtonian cosmologies to be spacetime manifolds {M., g), subject to the 
EFE, which contain a congruence of worldlines that is irrotational and shearfree. We call such a con- 
gruence a Newtonian-like timelike congruence, and denote the associated future-directed, normalised, 
tangent vector field by n (n a n a = — 1), the cosmologies by ( M., g, n). 

Most cosmologies will not admit such a congruence of Newtonian-like worldlines, and, hence, are not 
quasi-Newtonian; examples that are, are the FLRW cosmologies.^ If a cosmology does admit such a 
congruence, it will almost always be unique because it has a special relation to the Weyl curvature of 
( M., g ); we explore this further below. Because it is irrotational, this congruence is normal to a preferred 

FLRW geometry. 

4 And, in the non-cosmological context, the interior and exterior Schwarzschild spacetimes. 
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family of spacelike 3-surfaces of constant time, T: {t = const}, corresponding to the spacelike 3-surfaces 
of constant absolute time in a Newtonian spacetime. Because of the shearfree condition, these 3-surfaces 
are mapped conformally onto each other^Jby the preferred timelike curves n. The matter in such a setting 
will almost always not be comoving with these curves (or, equivalently, will not be at rest in the preferred 
spacelike 3-surfaces). The relativistic acceleration of these curves can be derived from a scalar potential 
that corresponds to the Newtonian gravitational scalar potential. 

In more detail: we aim to describe the physics of a pressure-free matter field from the perspective of 
a Newtonian-like class of Eulerian observers comoving with n. To this end, we apply a 1 + 3 covariant 
treatment adapted to n ( cf. Refs. [j?], Q |o) ) rather than the essentially equivalent ADM 3 + 1 
formulation of the dynamics of relativistic spacetime geometries [^[ The family of spacelike 3- 

surfaces T: {t = const} orthogonal to the timclikc reference congruence constitutes the rest 3-spaces 
of the Eulerian observers, with (time-dependent) 3-metric h ab := g ab + n a n b satisfying /i( a &) = = 
D a hbc- As to notation, for any geometrically defined variable T, we write T a ' b... := n c V c T a ■ and 
D a r b " c ... := h d a h b e . . . hs c . . . V^T 6 while angle brackets stand for a fully symmetric tracefree n- 
orthogonal projection as in, e.g., h^ ab ^ — [ h c ( a h d b ^ — i h ab h cd ] h c d- The timelike normals n to the T: 
{t = const} are irrotational: 

^ a (n)=0, (1) 

and shearfree: 

ffab(n) = . (2) 

The assumption (^|) corresponds to the vanishing of the anisotropic part of the extrinsic curvature of 
the orthogonal spacelike 3-surfaces T: {t = const} (these 3-surfaces exist because of property ([j]) ). The 
setting so outlined shares the features of the 'zero-shear hypersurfaces' gauge conditions introduced by 
Bardeen in his well-known work on linearised gauge- invariant cosmological perturbations |3| |] 

2.1.1 Zero magnetic Weyl curvature 

Through the 1 + 3 covariant Ricci identities for n JlJ, j2lfl , the kinematical restrictions given by Eqs. 
(|l|) and (|^) immediately lead to zero magnetic Weyl curvature of (M., g, n) as seen by the Eulerian 
observers comoving with n: 

H ab (n) = . (3) 

This has the important implication: there can exist no gravitational radiation in a quasi-Newtonian 
cosmology [p4[ . Being an invariant physical feature, it must also hold with respect to all other observers 
different from the Eulerian ones comoving with n. This shows that, considered as General Relativity 
spacetimes, such cosmological models are very special, for in fact generic motion of matter will generate 
gravitational radiation. However, its absence, of course, corresponds to the situation encountered in the 
Newtonian theory of gravitation. 



2.1.2 Acceleration scalar potential 

Preservation of the irrotationality condition (|l|) requires from the vorticity evolution equation (cf. Refs. 
Jill and ) that = e afcc Dftn c Q the spatial rotation of the relativistic acceleration n a (ri) :— n b V b n a 
of the timelike reference congruence n must vanish. Consequently, the acceleration is proportional to the 
spatial gradient of a 1 + 3 covariantly defined scalar field N jl5| : 

h a (n) := N- 1 B a N . (4) 

As already mentioned, this relativistic acceleration scalar potential corresponds to the Newtonian gravita- 
tional scalar potential. For convenience one may alternatively introduce a dimensionless scalar potential 
$ := ln(N/N ), where N is a constant.^ 

They are not mapped isometrically onto each other, primarily because we allow for an expansion of their normals. 
One could make a case that Strictly Newtonian- Like solutions should demand this stronger condition; however, that would 
exclude most of the expanding cosmological models that are a prime concern of this paper (but also are not strictly 
Newtonian). Accordingly, we stick to the weaker definition given above. 

6 The normals to the T: {t = const} are also non-shearing for Bardeen's 'longitudinal' gauge conditions fe|. 

7 In the comoving irrotational case of a barotropic perfect fluid, where the fluid acceleration is proportional to the spatial 
gradient of the total energy density, this is an identity. Here, however, it provides a genuine constraint on h a . 

8 The physical dimension of the lapse function N is [ length ] . 
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Taken together, the physical properties expressed by Eqs. (0) - (U) clearly justify the terminology 
'quasi-Newtonian' for the cosmological models at hand.f 

2.2 Local coordinates 

The relativistic acceleration scalar potential N is nothing but the prominent lapse function of the ADM 
3 + 1 formalism, which embodies (part of) one's personal choice of how to 'move forward in time' |53| 
within a dynamical spacetime setting. In terms of a set of (dimensionless) local coordinates {x" 1 }, the 
infinitesimal line element of a quasi-Newtonian cosmology can be written in the form (see, e.g., Ref. 

ds 2 = - N 2 (x l ) dt 2 + S 2 {x l ) da 2 , da 2 := / Q/3 (x 7 ) dx a dx (i , n a = N' 1 S a . (5) 

(the time coordinate is t = x°; spacetime coordinate indices i,j,k,... run through 0-3, and spatial 
coordinate indices ot,(3, 7, . . . run through 1 - 3), where the rate of expansion scalar O of the normals n 
is given by 

Q( X *) = 3-=3(SN)-'— , (6) 

and S, as usual, denotes a representative expansion length scale. Because this form of the line element 
implies that n is irrotational and shearfree, existence of local coordinates such that the line element has 
the form (^j is a necessary and sufficient condition that a cosmology is quasi- Newtonian in our sense 
(note that this form includes both the FLRW and Schwarzschild spacetimes). 

2.3 The matter flow 

The dust matter flow of our cosmological model^ is moving with average 4- velocity u (u a u a = — 1). 
Thus, relative to the Lagrangean reference frame comoving with u (where u a (u) := u b VbU a = 0), the 
energy-momentum-stress tensor takes the familiar form 

T ab (u) = flUaUb . (7) 

1 + 3-decomposing u with respect to the Eulerian n-frame we obtain |2^, |(| 

u a ■:= 7 (n a + v a ) , 7 := (1 - w 2 )- 1 / 2 , v 2 :=v a v a >0, v a n a = 0. (8) 

The variable v a is the peculiar velocity of the matter flow with respect to the T: {t = const}. Hence, the 
1 + 3 covariant decomposition relative to the Eulerian n-frame ( where U a b := —n a nb, h a b = 5 a b — U a b ) 
leads to 

T a b(n) = l*>n a n b + 2 q {a n b) +ph ab + ir ab ; (9) 
the matter variables are given by p^ j p 2 ] 

[i = ^ 2 Ji, p=±nv 2 , q a = fiv a , ir ab = (1, «(„ v b ) . (10) 

Note that p and n a b, which are typically non-Newtonian sources of spacetime curvature, are quadratic 
in v a , so they will be negligible when we linearise our models about a non-comoving FLRW background, 
assuming that v a itself is a deviation variable of first-order smallness. However, q a is linear in v a ; this 
will be crucial later on. The dependence of the fluid kincmatical variables a ab , 6 and Co a b on the set 
{n a , n a , 9, v a , D a Vb } is given explicitly in subsection A. 3 of the appendix.^] In the limit v a — > the 



situation at hand reduces to the standard comoving description of FLRW cosmological models with dust 
matter source field (in this limit, the acceleration of the normals has to vanish, because the dust flow 
lines are geodesic). 

9 For further discussion sec Ref. and the discussion of 'Newtonian-like' spacetimes in Ref. |3. 

10 In the case when S = S(t), we can alternatively use the conformal time coordinate rt(t) = I dy/S(y) of physical 

dimension [length] -1 , integrated along the integral curves of the timelike reference congruence n, and the line element 
takes the form ds 2 = S 2 (r])( — N 2 (x l ) drj 2 + da 2 ). Note, however, this does not work in general when S = S(t, x a ); for then 
r\ = rj(t, x a ) and there will be time— space cross terms in the resultant form of the line element. 
In our discussion, for simplicity, we set the cosmological term to zero, A = 0. 

12 We could consider several matter components: each would obey similar equations, with the total matter variables simply 
being the sum of those for each component. Note the direct relation to the approach taken in formulating thermodynamical 
properties of matter within kinetic theory. 

13 For the covariant derivatives of n a and v a one obtains, respectively, V a n(, = — n a rtj, + g©/i a (, and V a li(, = 
- (v c h c ) n a n b - n a V(fi) + | Ov a n b + D a v b . 
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2.4 Uniqueness of N and n 

As is well-known, the ADM 3 + 1 formalism does not provide one with a natural evolution equation for 
the lapse function N (nor the shift vector N a ), and, therefore, by Eq. (^), also not for the acceleration 
h a of the normal congruence n f45|| . Hence, it is a major dynamical issue as to how to determine, 
for a given physical problem, the best prescription for moving forward in coordinate time t from one 
spacelike 3-surface within the foliation to the next In general, a 'time slicing condition' leading 

to a typically differential constraint relation has to be introduced by handj^] However, a fundamental 
point has already been raised above: despite this lack of uniqueness of N from a dynamical viewpoint, 
we expect the normal congruence n in our case to be highly restricted geometrically , indeed to usually 
be uniquely determined, because of the conditions (a) that Eulerian observers comoving with n will see 
zero magnetic Weyl curvature, Eq. (||), and (b) have a highly restricted electric Weyl curvature, or, 
equivalently, a very restricted intrinsic curvature allowed for the spacelike 3-surfaces in the family T: 
{t = const} (this will be explored below). The hope then is that choice of this congruence essentially 
determines the time evolution of N, because it determines the acceleration h a ( which occurs in Eq. (jij) ); 
so uniqueness of n will imply a unique choice of the lapse function. 

Consequently, it is important to determine when the Newtonian-like timelike reference congruence n 
really is unique. As regards point (a) above, if we express the Weyl curvature relative to n in terms of 
its electric part, we obtain 

C ah cd (n)=^ a n [c E b \+Ah^ [c E^ d] . (11) 

Hence, the magnetic Weyl curvature relative to another timelike congruence, n, where n a = 7 (n a + v a ) 
(cf. Eq. (|)), is 

H ab {h) = 2 7 2 e cd{a [ n b) {E c e v e ) + E c b) ] v d . (12) 

So, generically, for E ab + 1 and v a 7^ 0, we have H ab ^ 0. On the other hand, both H a b{ n ) — and 
H ab (n) = 0, iff 

tcd(a [ n b) (E c e v e ) + E c b) ]v d = 0. (13) 
This will be satisfied in the FLRW case, where = E ab (n) = H ab (n) (implying = E ab (n) = H ab (n); cf. 



subsection A. 3 of the appendix), and also for cosmologies with even shearing timelike reference congru- 
ences that possess LRS symmetry if v a characterises a Lorentz-boost along the 1 + 3 covariantly defined 
preferred spacelike direction e Jll], ^2j. Equation ( |l3| ) is a necessary (but not sufficient) condition that 
the same spacetime geometry can appear Newtonian-like relative to two (or more) timelike congruences. 
When this is satisfied, one then also needs to find if (b) is satisfied: does the electric Weyl curvature 
have the highly restricted form that allows a zero initial rate of shear to remain zero, given the matter 
source fields present? These conditions are rather complex, nevertheless, we would like to get a complete 
solution: 

Problem 1: Determine all quasi-Newtonian cosmologies that allow more than one 
Newtonian-like timelike congruence. 

This family includes constant curvature spacetimes and the FLRW spacetimes; are there any others? 



3 1 + 3 covariant dynamical equations 

Insisting on the timelike reference congruence n to have zero rate of shear, o~ ab = 0, converts the original 
a-( ab ) -equation (a Ricci identity determining the evolution of a ab along n) into a constraint equation which 
we can regard as either restricting the allowed time slicing, given the Weyl curvature of ( AA, g, n) (i.e., 
as a differential time slicing condition for the lapse function N), or as restricting the allowed electric 
Weyl curvature, E a b(n), given a choice of the time slicing of (M, g, n). In either case, the constraint so 
generated assumes the explicit formp| 

= E ab - D< V> - n< a n b > - \ fj, v {a v b) . (14) 

14 This indeterminacy does not arise when the congruence of preferred worldlines is comoving with well-specified matter 
(e.g., a perfect fluid with specified equation of state). However, that is not the case here. 

15 We draw the attention of the reader to the fact that the differential time slicing condition constituted by Eq. (M) is con- 
ceptually analogous to obtaining a relativistic generalisation of the PFE when imposing 'maximal' time slicing characterised 
by© = 0=>6 = (but c a (, ^ 0) in the case of asymptotically flat spacetimes Wq] . 
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This is a highly restrictive assumption; we wish to explore its implications. 

Using the variables specified in Eq. ([H]), under the stated assumptions ([!]) - characterising 
quasi-Newtonian cosmologies with dust matter source field one derives from the general 1 + 3 covariant 
dynamical equations (cf. Refs. Q and [^l] ) the set: 

Evolution equations: 

6 = -i9 2 + D a n a + (n a n a )-i(l + trV (15) 

fi = -(l + ±v 2 )efi-v a D a fi~LiD a v a -2Li{v a n a ) (16) 

u<°> = -i(l-i; 2 )e/ -v b D b v a + (v b n b )v a - h a (17) 

+ | (D c i; c ) jj, v {a v b) - \n D<V> + n v c v {a D ( V> + n v c v {a D [ V ] . (18) 

The evolution equation for v a is a consequence of the original g^-equation (Bianchi identity); it estab- 
lishes conservation of momentum (density) . As such it is the relativistic analogue of the Euler equation 
in the Newtonian theory. 

Primary constraint equations: 

= (Ci) a := D a e-^v a (19) 

= (C 2 ) Q := e abc B b n c (20) 

= (C 3 ) ab := £ Qb - D <a n fc > - n< Q n b > - §/iu< a u 6 > (21) 
= (C 4 ) a := D^" 6 - I (1 + I w 2 ) D> + i [ 9 + | /i- 1 (/D 4/1 ) + | (D b v b ) }[iv a 

+ ip 6 [D ( V>-5D [ % t l]. (22) 

The constraint (Ci) a , being the reduced form of the original D{,<7 ab -constraint (Ricci identity), states 
that the peculiar motion of the matter as quantified by v a induces a non-zero spatial gradient in the 
expansion rate of the normals n. Constraints (C^)" and {C 3 ) ab are the key new constraints resulting 
from the quasi-Newtonian assumptions, (C^)" leading to the existence of the relativistic acceleration 
scalar potential N (see Eq. (Q) ). {Ci) a is the D;,i? ab -constraint (Bianchi identity). Additionally, there 
are constraints related to H ab : 

= (C H ) ab := H ab (23) 
= (CW) a := D b (C H ) ab + ±e abc D b ( f iv c ) . (24) 

(CH) ab holds as long as Eqs. (||) and (||) apply. Using Eq. ( g6| ) of subsection Al in the appendix, 
the (Co.tf ) a -constraint (Bianchi identity) can be reduced to = D b (CH) ab ~ \ e abc T) b {Ci) c and thus is 
identically satisfied if each of (Ci) a and (C#) ab is valid. It then follows that 

e abc T) b v c = e abc v b ^ & D [a v b] = ^ v [a D b]l i ; (25) 



the spatial rotation of v a is equal to the algebraic cross product between v a itself and the fractional total 
energy density gradient (D a /i)//i. As such, the vorticity of v a will be of second-order smallness when we 
only consider first-order deviations about a FLRW background setting associated with the T: {t — const}. 
Altogether, the fully orthogonally projected covariant derivative of v a can now be decomposed as 

D a v b = D <a w b ) + | (D c v c ) h ab + /.r 1 v [a TJ b] fi ; (26) 

D(ai'fe) and D a w a representing the distortion (or shear) and spatial divergence (or expansion) of the 
peculiar velocity field, respectively. With condition ( p3| ) prescribed, the status of the -equation 
(Bianchi identity) can take either the rank of yet another constraint relation (cf. Ref. p4|] ), or it 
naturally yields the propagation of (Cy) ab along n. Here we find that 

{C k ) ab = (Cff)< o6 > = - 9 (C H ) ah - [ D <a + n <a ] (C 2 f - e cd ^ a [ D c + n c } (C 3 )% ; (27) 
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if {C2Y and (C^)" 6 are satisfied on %: {to = const}, then a consistent set of initial data for quasi- 
Newtonian cosmologies with dust matter source field characterised by Eqs. ([j]) - (Q) generates no magnetic 
Weyl curvature along n. 

A particularly important restriction on (A4, g, n) results as a compatibility requirement between 
(C 3 ) ab and (Ci) a . It follows that 

= D b (C 3 ) ab = (C 4 ) a - n b (C 3 ) ab - § (C 5 ) a , (28) 

where (Cs) a is given by 

(C 5 ) a := D a [D b n b + {h b n b ) } + n a [V b n b + {n b ii b ) ~ \Q 2 + {\ - \v 2 ) ^] 

- i (1 + v 2 ) D> + i [ 9 + 3 ^ (v b B bt i) + f (D b v b ) + 3 (v b h b ) } p v a 

+ ^v b [B< a v b) -5D^ a v b] } , (29) 

and = (Cs) a is demanded. Using Eq. (|J), this provides a third-order partial differential equation for 
N. 

In summary, the particular non-trivial constraints characterising the quasi-Newtonian case at hand 
are (C 3 ) ab and its immediate consequence (C5) , for these embody the consequences of the zero-shear 
requirement (|^). 

3.1 The geometry of the orthogonal spacelike 3-surfaces 

Next, we give the geometrical variables characterising the intrinsic curvature properties of the orthogonal 
spacelike 3-surfaces T: {t = const}, which exist because of the requirement ([!]). 

Trace and tracefree parts of Gaufl embedding equation: 

= (C G ) := 3 i?+§9 2 -2/i (30) 
= {C G ) ab := 3 S ab - E ab - \ fjLV {a v h ^ . (31) 

3 R and 3 S a b denote the respective irreducible parts of the 3-Ricci curvature. As demonstrated for the 
example of a Lagrangean description of an irrotational dust matter source field in Rcf. , the GauB 
embedding equation has to be included in order to obtain a set of (1 + 3 covariant) dynamical equations 
which is equivalent to the EFE. Note that the geometrical structure of Eq. ( p0| ) is identical to the 
Friedmann equation in FLRW cosmology. 

3- Cotton- York tensor Q: 

3 C ab := h 1/3 e cd{a D c 3 S b) d 

= -2^/V d < Q [ n c D^n d) + ± M « c (T>^v d) - D^v d] ) } - h 1 ' 3 ^ v< a e^ cd h c v d 

+ h 1 ' 3 e cd < a [ B c (C G f d + B C (C 3 )\ - n c (C 3 f d ] + h 1 ' 3 h^ a (C 2 f . (32) 

The 3-Cotton-York tensor, although not having the dimension of a curvature variable, encodes the 
conformal curvature properties of the T: {t — const}. They are conformally flat if 'KJab is zero. It has 
been suggested that a non-zero 3-Cotton-York tensor could be intimately related to the presence of 
gravitational radiation in a spacetime geometry (see, e.g., Ref. p2[ ). Consequently, one would expect it 
to vanish in the quasi-Newtonian case. 

3.2 Propagation of and h a 

Having now attained the full set of 1 + 3 covariant dynamical equations for quasi-Newtonian cosmlogies 
with dust matter source field, it is clear (in accordance with the ADM 3 + 1 results) that there is no 
explicit equation determining the evolution of N or n a . Two differential expressions containing the latter 
particularly spring to the eye. In terms of N the ri a -part of Eq. (^) is given by 



D<a«b) + n {a h b) = N 1 ~D {a D b) N 



(33) 
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The symmetric tracefree derivative D^h^ is the distortion of h a . On the other hand, the spatial diver- 
gence T) a n a can be interpreted as its expansion. Again, in terms of iV we have 

D a ri a + (n a n Q ) = A" 1 D a D a A ; (34) 

this is a source term in Eq. ( |l5| ) and corresponds to the Laplacian of the gravitational scalar potential in 
the Newtonian theory. We also have no evolution equations for neither of these two differential quantities. 
This remains true if one writes out the full set of equations in an associated 1 + 3 orthonormal frame 



formalism, which we give in subsection kA of the appendix. Hence, we face 



Problem 2: In those cases where the Newtonian-like timelike congruence is unique, how 
is the corresponding evolution of the lapse function N along these worldlines determined? 

Lacking explicit evolution equations for any of the relativistic acceleration^ its spatial divergence p| 
or other spatial derivatives, we presume the former must be implicit in the consistency conditions that 
arise when we impose the quasi-Newtonian restrictions (|l|) and @ on our timelike reference congruence. 
However, we have been unable to obtain a general prescription as to how this happens. 

3.2.1 An Ansatz 

In a formal sense, an evolution equation for n a in the quasi-Newtonian dust case can be derived by 



combining Eq. (^) with the commutator relation (|92|) in subsection A.l of the appendix, thus yielding 



(n)'<"> = N~ 1 D a (N) - \&h a . (35) 

The latter makes explicit the necessity for an evolution equation for N, required to close the set of 
dynamical equations. Failing any other way of determining this, we propose to evolve N along n according 
to the Ansatz 

N = a 9 N where a — const , (36) 

which is a common choice when modelling dynamical spacetime geometries by numerical methods ; for 
a = 1, e.g., this Ansatz specialises to the so-called 'harmonic gauge' where N propagates along n via a 
wave equation at the speed of light (cf. Refs. and [|| ). Other particular values for the parameter a 
will be mentioned in section || below. With this assumption, it follows that 

(n)' <Q> = (a - l)Qn a + a^iv a + a(d) a . (37) 

Thus, the evolution of the dynamical quasi-Newtonian cosmologies is now deterministic. 



4 Constraint analysis 

To check whether quasi-Newtonian cosmologies with dust matter source field can provide consistent 
solutions of the EFE requires a detailed analysis of the preservation properties of the constraint relations 
as a model evolves. That all constraints involved be preserved along the integral curves of n is a necessary 
condition for consistency of the 1 + 3 covariant dynamical equations; a full consistency check involves 
additionally showing, e.g., that all Jacobi and Ricci equations of an associated 1 + 3 orthonormal frame 
formalism are satisfied ( see, e.g., Refs. jllj and j2^| ), or casting the evolution equations into a first-order 
symmetric hyperbolic form (see, e.g., Refs. |j| and for the cases of a vacuum and a perfect fluid 
spacetime geometry, respectively ) . 

16 This contrasts the case of a barotropic perfect fluid represented in a comoving frame, where (i) the relativistic acceler- 
ation is determined from the spatial gradient of the pressure by the momentum conservation equation, (ii) the pressure is 
determined from the total energy density via the equation of state, and (iii) the energy density time evolution is determined 
by the energy conservation equation. In the present situation, because the reference frame is not comoving, step (i) fails; 
given any specified acceleration distribution, the momentum conservation equation determines the peculiar velocity of the 
matter fluid relative to the Eulerian reference frame. 

17 In the strictly Newtonian case we would additionally have © = 0, and then the Raychaudhuri equation ( p^ ) becomes 
the relativistic generalisation of the PFE determining this spatial divergence. In the case of an expanding cosmological 
model, however, this equation rather determines the time rate of change of the expansion that is compatible with whatever 
value holds for this spatial divergence. 
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As mentioned, ideally we would like to derive the form of N in the generic case from this consistency 
analysis, but have been unable to do so; consequently, a prescription for N is needed, followed by a 
subsequent consistency check. We assume the evolutionary behaviour of the lapse function N along 
the integral curves of n to be given by the Ansatz (J3q ) . Using relations for the propagation of spatial 
derivative terms given in Ref. PU P and also Eq. (|96|) in subsection A.l of the appendix, we then find 
that the primary constraints fll9|) - (22|) and (29) propagate according to 

(Cx)^ = -eidr+iar (38) 

(C 2 )^ = ( a -|)6(C 2 ) a + ie Qbc n b (C 1 ) c (39) 
(C 3 ) {ab) = -e(C 3 ) afc -aD< Q (Ci) b> -4(a- ^)n< Q (Ci) b> 

- (a + |) [ 9 D ( V> + | M D ( V> + | u<° D b) M ] " 3 (a - ±) 6 n< a n b > 

- 6(a- ±)/W a 7j b > +2(id c i>< 11 [ D^/* + D^V ] 

+ [ (1 - if 2 ) 9 + /i- 1 (« c D c /i) + f (D c t/) ] M u< a v b > (40) 

(&) (0) = -|e(c 4 ) a + i(i + « 2 )M(Ci) a -^f a [^(Ci) b ]-i? Q fc (c 1 ) b + i/i(c 6 ) Q (4i) 

(C 5 )< a > = (a - 1) 9 (C 5 ) a + a [ D a + 2 n a ] [D b (d) b ] + 4 (a - ±) [ D a + 2 h a ] [n fc (d) b ] 
+ 2 [ (a - |) TJ b h b + 2(a-±) (n b n b ) - a \ 9 2 

+ i( a + i( tt -2> 2 -y)H(Cir 

+ !(«-§ + 5« 2 ) K(Ci) b ] + § (a + | + lv 2 )^D a {D b v b ) 
+ \ {a + i + (a + f> 2 - §u 4 ) 6 D> - (a + |) S 2 ^v a 

+ | (a + | + (a - ^)w 2 + |w 4 ) /i 2 u a + further non-linear terms , (42) 

where {C§) a is defined by 

(C 6 ) a := D Q (D b 7j b ) - §D b D<V> 

- [ i 6 2 - (1 - ft; 2 + f§« 4 ) M + f 9 (D b t, b ) + 2 « b D b (D cU c ) + | fj,- 1 (D b D b M ) 

+ 3 (u^DM) - 3 (D [6 u c] DM) + 3 (v b h c D [b v c] ) 

+ 3^ _1 (wiDMD c /i) ]v Q 
+ w b [ 3 7j c D c D [ % b l + | /i- 1 D b D a /i + 5 (D c t; c ) D [ VJ + 3 fx' 1 (w c D c /i) D [ VJ 

+ 3 (v c n c ) D [Q 7j b l + | D< V> + § n< a n b > ] 
+ \ (D b w b ) /i- 1 D a /i - | D< V> /i- 1 (D 6/ x) - \ DM /i- 1 (D b/ i) 
+ 3 D< V> D [b 7j c] v c + 3 D [ V ] D (b 7j c) t/ - 6 D [ V ] D [b 7j c] v c . (43) 

We remark that, besides {C\) a to (Cs)", in principle also (Cg) and (Cc) ab given in Eqs. © and © 
are spatial constraints which need to be preserved along n. However, the 1 + 3 covariant formulation 
does not provide natural evolution equations for the spatial curvature variables 3 R and 3 S ab , or, rather, 
the underlying connection components relating to h ab . To complete the constraint analysis, one instead 
needs to introduce either a 1 + 3 o rthon ormal frame ( and the equations deriving from the Jacobi and 
Ricci identities listed in subsection A. 4 of the appendix), or a set of local coordinates (see, e.g., Ref. 
ijllf ). Experience shows, nevertheless, that in general both (Cq) and (Cc) ab are preserved (cf., e.g., 
Rcfs. H and |l| ). 

Returning to the set © - ©), we now see that with the choice © both (Ci) a and (C 2 ) a will be 
preserved along n, if each of (Ci) a , (C^)" and (Cs) a hold on an initial spacelike 3-surfaces Tq: {to = const}. 
However, propagation of (Ci) a generates the secondary constraint (C6) a .0 Furthermore, it is clear that 
neither (C?,) ab nor (C$) a can be made to be preserved along n for a single value of a, so they also 
generate secondary constraints. Hence, exact quasi-Newtonian cosmologies with dust matter source field 

18 Note the error in Eq. (B.13) of Ref. ]2l|: the numerical factor preceding the last term on the RHS should be '3' rather 
than '2'. 

19 Note that the structure of (Ce) a is quite reminiscent of the general D^o^'-constraint (Ricci identity) in that it relates 
the spatial divergence of the shear D^v^ of the peculiar velocity field v a to the spatial gradient of its expansion D a v a . 
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and a lapse function evolution given by Eq. ( pq ) are not consistent at the first derivative level as regards 
propagation of the constraints along n — new conditions arise. The question is whether consistency can 
be attained at a higher derivative level, or with a prescription for N alternative to Eq. (|36|); however, 
in view of the complexity which arises in the situation just discussed this seems doubtful except in very 
special cases. Indeed, initially our hope was that this project could lead to a substantial class of exact 
non-linear general relativistic generalisations of the Zel'dovich approximation method for gravitational 
instability scenarios; however, considering the complexity of the consistency equations, and noting the 
physical restrictions implied by excluding all gravitational radiation phenomena, we now believe there 
exist only very few consistent quasi-Newtonian cosmologies with dust matter source field. 

We have been unable to resolve the generic consistency question by either 1 + 3 covariant or orthonor- 
mal frame methods. Thus, we end up with the open issue: 

Problem 3: Determine (a) the set of cases where the constraint equations for a quasi- 
Newtonian cosmologies with dust matter source field are consistent with the time evolution 
equations, and (b) the related set of consistent quasi-Newtonian cosmologies. 

In the following section we derive some of the consequences that arise from imposing particular dynamical 
restrictions on quasi-Newtonian cosmologies with dust matter source field. Thereafter, in section |[ we 
linearise the dynamical equations about a non-comoving FLRW background. 



5 Constraint analysis in restricted cases 

We have examined the reduced set of 1 + 3 covariant dynamical equations when particular geometrical 
restrictions are placed on the quasi-Newtonian cosmologies with dust matter source field, to see if we can 
solve the problem in those cases. In particular, we have examined quasi-Newtonian cosmologies with (a) 
zero spatial rotation of v a , (b) zero spatial divergence of v a , (c) LRS spacetime symmetry. In each case 
we find a more manageable reduced set of equations, but nevertheless have been unable to resolve the 
problem. 

Where we have made more progress is in the cases of quasi- Newtonian cosmologies with (a) zero electric 
Weyl curvature (the conformally flat case), (b) the comoving case, and (c) zero relativistic acceleration. 
We now consider those cases in turn. 

5.1 Zero E a b — FLRW in a non-comoving frame 

Specialising to zero electric Weyl curvature with respect to n, it follows from Eq. (|ll|) that 

= E ab = H ab o C abcd = . (44) 

Now for fj, ^ and 9 ^ it is well known that the only conformally flat non-static dust cosmologies belong 
to the FLRW family ( see, e.g., Refs. |7, 12, 29[ [35| ). Thus, the expanding non-empty cosmological models 



of this kind are just the FLRW geometries as seen by a non-comoving Eulerian observer in a Newtonian- 
like rest frame (in this case these frames are not unique, because of the vanishing Weyl curvature, cf. the 
discussions in subsections 2A and |6.2.4[ ). 



From Eq. ( J22| ) we then find that (64)" becomes 
(C* 4 ) a = -| (1 + \v 2 ) D a fi + i [ 9 + § m" 1 (v b D bM ) + f (D b v b ) }fiv a + ^v b [ D<V> - 5 DM ] , (45) 

(essentially: the equation D b E ab = 0), which can be treated either as an equation for D a /i (determining 
the change of fx along the direction v a in T: {t = const}),0 or for T) a v b . Then, Eq. ( |l"g| ) provides the 
new constraint (essentially: the equation E^ ab ' = 0) 

{C) ab := § (1 - it- 2 ) 9 ^ + i (/D c p) t)<" w fc > - i «< a D 6 V 

+ l(D c v c )[iv {a v b) - \nV {a v b) + f iv c v {a ~D {b) v c) +nv c v {a TJ [b) v c] , (46) 

which is also an equation for D a /z or D a tv We will discuss the non-comoving FLRW case in linearised 
form in detail in section ^| below. 

20 Even though this is a FLRW model, this will be non-zero if v a ^ because then we are in a tilted frame ( cf. Ref. 
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5.2 Zero v a — The comoving case 

If we impose the comoving condition v a = 0, it immediately follows from Eq. ( p~7[ ) that the relativistic 
acceleration vanishes P] 

v a = n Q = . (47) 

Then the equations simplify dramatically because by Eq. (|2l| ) 

^ab = . (48) 

Thus, as condition ( fli]) is satisfied, the only cosmologies satisfying Eq. ( [47j ) are the FLRW models in 
their standard (comoving) frame. 

In detail: as n a — 0, we have by Eq. (||) that N = N(t); the relativistic acceleration scalar potential 
is spatially constant Q and we can choose N = 1 without loss of generality such that t in Eq. (||) then 
measures proper time of physical dimension [length] (see, e.g., Ref. [|l2| ). From (Ci) a = 0, D a 9 = ^> 
we can choose local coordinates so that S — S(t) in Eq. (||). From (62)™ = 0, D a /i = 0, so Eq. ( |30| ) 
=4> D a 3 i? = and Eq. ( |3l| ) ^afc = 0. Hence, these cosmologies are spatially homogeneous (i.e., 
/i = n(t), O = O(i), 3 i? = 3 R(t)) with orthogonal spacelike 3-surfaces with 3-metric h a b of constant 
curvature K(t) := 3 i?/6. Thus, the 3-metric fapix 1 ) in Eq. (||) can be chosen in one of the standard 
forms for a 3-space of constant dimensionless scalar curvature k where k can be normalised to one of 
{ 0, ± 1 } by suitably renormalising S(t) ( see, e.g., Ref. ]14| ). The primary constraint equations are then 



all identically satisfied, and the remaining non-trivial propagation equations are Eqs. ( |15| ) and (16), with 
a first integral given by Eq. (|30|), provided K(t) = k/S 2 (t), which will indeed hold. 

It is fundamentally important for what follows later, that by condition (C\) a = the comoving case 
is precisely that case where the spatial gradient of the expansion rate is zero: 

v a = <S> q a = & D a 9 = . (49) 

Thus, we note that any of these conditions imply a comoving FLRW subcase of our quasi-Newtonian 
setting. 

5.3 Zero h a — The geodesic case 

When the relativistic acceleration of the shearfree normals n vanishes, again 

ii a = ^ N = N(t) ; (50) 
the preferred (Newtonian-like) Eulerian observers are in free fall, and so experience no effective gravita- 



tional field. The canonical choice for the lapse function is N = 1 (cf. subsection ^2). Indeed, we have a 
zero-acceleration quasi-Newtonian cosmologies if and only if there are local coordinates in which the line 
element reduces to Eq. ([|) with N = 1. Examples are Minkowski spacetime and the family of FLRW 
models. The question is whether there are other quasi- Newtonian cosmologies satisfying Eq. (^p|). For 
the conditions specified by the latter, Eq. (|35|), which does not assume a particular form for N, only 
leads to an identity and so provides no new information. On the basis of the Ansatz ([36]), however, we 
find from Eq. @ that for a ^ 0, /i ^ and (d) a satisfied, 

h a = v a = ; (51) 

this is just the comoving case just discussed in subsection eT3. 



6 Linearised perturbations of a FLRW dust cosmology 

When discussing the evolution of linear spatially inhomogeneous deviations from a FLRW spacetime 
geometry, it is a common feature of many papers to introduce local coordinates and an explicit Ansatz 
for the infinitesimal line element (see, e.g., Refs. |j, ^, [36| ). However, in doing so important physical 

21 We have imposed the condition that the matter is pressure-free (and so moves geodesically), and now that matter is 
moving along the integral curves of n; hence, these must be geodesic. 

22 That is, it is constant in the preferred spacelike 3-surfaces T: {t = const} orthogonal to n. 
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aspects of the dynamics of the setting often become rather obscured. We claim that it makes for clarity 
to work in 1 + 3 covariant terms for as long as possible. In this section we discuss the important 
issue of a consistent Newtonian limit of FLRW-linearised cosmologies, once we have provided a 1 + 3 
covariant formulation of the setting at hand. For a comparison discussion, we first briefly review a recent 
treatment given by Bertschinger of the evolution of linear energy density inhomogeneities (i.e., so-called 
'scalar modes') as induced by matter in non-relativistic motion Q. The reason for referring to his work is 
that the timelike reference congruence he chooses is of the Newtonian-like kind, i.e., it is irrotational and 
shearfree. Other well-known discussions such as, e.g., that given by Peebles in Ref. |36| pick the timelike 
reference congruence to be irrotational and geodesic instead ( in Bardeen's terminology, the latter choice 
provides 'synchronous' gauge conditions § ). We briefly refer to them at the end. 



6.1 Bertschinger's method 



Bertschinger chooses what he calls 'conformal Newtonian' gauge conditions Q. In Bardeen's terminol- 
ogy these are 'longitudinal' gauge conditions, where the timelike normals n to the spacelike 3-surfaces 
T: {t — const} are shearfree || , and it is assumed that the matter perturbations do not produce any 
anisotropic pressures. The discussion in Ref. Q] then focuses on scalar linear spatially inhomogeneous 
deviations. Choosing a set of Cartesian spatial coordinates x := (x,y,z) comoving with n on a spatially 
flat FLRW background as well as a conformally rescaled time variable 77 defined by^ 



d 
dt 



1 d 
S dn~ 



dt 2 



1 

S 2 



d 2 1 dS d 
dr\ 2 S dr\ dr\ 



(52) 



it is stated that only one (scalar) variable, <&, is needed to express the deviations in the background line 
element. Though no explicit explanation is given, this follows from the assumption of zero anisotropic 
pressure perturbations in combination with part of the linearised EFE ]3^| ; we give the relevant details 
in subsection [A.5| of the appendix. One so obtains the line element reduced to the simple form 



ds l 



S 2 (n) [ - (1 + 2$) drj 2 + (1 - 2$) {dx 2 + dy 2 + dz 2 ) 



(53) 



No further statement is made in Ref. ||] as to the 77-dependence of <&. 

In the local coordinates employed, the proper peculiar velocity of a matter particle is given by v = 
dx/drj. It should here be stressed that in the Eulerian description adopted by Bertschinger the peculiar 
velocity "is defined as the average momentum per unit mass of the particles in the vicinity of a given 
Eulerian position" Q, i.e., in terms of the momentum density of the non-relativistically moving matter. 
Gravitational effects generated by the momentum density, however, are said to be neglected. Dynamical 
equations are now derived from the remaining linearised EFE on the basis of the assumptions that (1) the 
observable part of the Universe is nearly FLRW on scales of the present day Hubble radius, (2) | $ | <C 1 
and v 2 <C 1, (3) only spatial inhomogeneities on scales much less than the present day Hubble radius 
are considered, and (4) the deviations in the energy density are dominated by non-relativistically moving 
matter only, so no anisotropic pressures need to be accounted for [Q. For dust, the set presented in Ref. 
(I to evolve and constrain the deviations corresponds to 



<y = 







S' 
3-^ 



fi m V • v - V • (5/j, v) 



- — v- (v ■ V) v — V$ 
V • V$ - \ S 2 6(i ; 



(54) 

(55) 
(56) 



the prime denotes a partial derivative with respect to conformal time 77, the V-operator acts in terms of 
the comoving Cartesian coordinates x, and the definition S/u, := fj, — fi m was used (jj, m being the mean 
total energy density of the non-relativistic matter) . Equation (^4j) is the continuity equation representing 
conservation of energy (density) , Eq. ( |55| ) is the Euler equation representing conservation of momentum 
(density), while Eq. ( p36| ) is the PFE of Newtonian gravitation. In view of the Ansatz ([53|) for a spatially 
flat background, the PFE describes how matter inhomogeneities induce the irregularities in the almost- 
FLRW line element. 



'If, contrary to the time coordinate in Eq. feh, t here denotes dimensional proper time, then 77 is dimensionless. 
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Keeping only terms of first-order smallness, Bertschinger derives from the set ( p4|) - (j56|) the standard 
second-order linear ODE [Q, |3(| 

= 6" + ^S' -\S 2 [i m 6 (57) 

for the dimensionless density contrast S := Sfi/ /i m . This is, however, a highly gauge and coordinate depen- 
dent procedure. By contrast, a corresponding equation can be derived in a more gauge-independent and 
covariant way by using geometrical variables and their dynamical equations. In the comoving case, this 
reproduces exactly the equation (|57]) (see Ref. fl7|| ); however, in the non-comoving case corresponding 
to the situation discussed here, generically it gives a different result, as we now explain. 



6.2 Covariant and gauge-invariant formulation 

We contrast the description given by Bertschinger with a non-comoving 1+3 covariant and gauge- invariant 
formulation of linearised perturbations of FLRW models in the spirit of Ellis and Bruni 1 17 , but in the 
quasi-Newtonian cosmology context. That is, (consistent with Eq. ([53|) ) we assume the time slicing 
of (M, g, n) provided by the spacelike 3-surfaces T: {t = const} with shearfree normals n to have an 
approximately FLRW geometry, such that the quantities v a and h a can be regarded as (gauge-invariant) 
deviations of first-order smallness. As we consider non-relativistically moving matter only, we do account 
for its momentum density q a , which is linear in u a , but neglect its isotropic and anisotropic pressures, p 
and 7r a fc, as they are of second-order smallness in v a ( cf. Eqs. ( fl0| ) ). The difference to the Ellis-Bruni 
treatment is that here we are not using a comoving choice of the timelike reference congruence, thus, 
the procedure is frame-independent iff the Newtonian-like timelike reference congruence n is unique. We 
expect this to be the case in general, but, as discussed in subsection |3l| above, there may be some freedom 
in its choice in specific restricted geometrical cases. In particular, we can view a strictly FLRW model 
from a tilted (non-comoving) reference frame. We will come back to the aspect of factoring out the pure 
non-comoving FLRW subcase, once we have investigated the general properties of the FLRW-linearised 
quasi-Newtonian dynamical equations. 



6.2.1 The lapse function in almost-FLRW models 

Given this context, we now return to the question as to how to evolve the lapse function N along the 
integral curves of our timelike reference congruence n. The standard literature often deals with this 
aspect by merely stating (though generally not providing explicit mathematical relations; however, see, 
e.g., Refs. 0, 1?], H for cases where the background is expanding, and Ref. [[l4| where it is not) that 
the non-relativistic peculiar motion of the matter induces only 'slowly' evolving variations in the metric 
structure of ( A4, g, n). One way to interpret statements of this kind is to take N = 0, but other options 
may also be taken into account. A number of possible ways of evolving the lapse function N along n 
are contained in the Ansatz ([56|), which we continue to use in the subsequent analysis, leading to the 
evolution equation ([3?]) for h a . Besides the 'harmonic gauge' given by a = 1, this Ansatz also allows for 
the special subcases: 



(i) a = — j (SN)' = 0; the dimensionless comoving (with the Eulerian observers) lapse function 
(SN) is covariantly constant along n; this choice is suggested from computing to first order the 1 + 3 
ONF commutation functions associated with the line element ([53]): as can be derived from the relations 

= o a fi = u a = n a p 



given in subsection A. 5 of the appendix o ne o btains h a 



f2 a , and so the 



evolution equation for n a is given by Eq. ( 124 ) in subsection A. 4 



(ii) a = => N = 0; as mentioned above, here the lapse function N itself is covariantly constant along 



(iii) a = 4; provided one makes the identification N = 5(1 + 2<&) 1 / 2 J^| this value is implied by 
the analysis of growing mode 'adiabatic' perturbations in spatially flat FLRW backgrounds given by 



24 At this point, one faces a problem with the physical dimension of the time coordinates involved. While our t in Eq. (g) 
is dimensionless, letting (as is standard in the ADM 3+1 formalism) N carry the dimension of the line element, in other 
works such as Refs. E| and p?| t is dimensional proper time (and N = 1). This prevents an unambiguous identification 
of N with variables such as S and <I>, once a (dimensionless) conformal time coordinate 77 is introduced. 
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Mukhanov et al Q, p225, where (in their notation) $' = 0. 
6.2.2 Evolution of constraints 

Besides v a and n a , also the 1 + 3 covariant quantities X a := D a /x and Z a := D a 9 encode deviations of 
first-order smallness from an exact FLRW geometry of our time slicing; indeed we will regard X a and Z a 
as our prime perturbation variables. The evolution equations for these variables are 

FLRW -linearised •perturbation evolution equations: 

X^ = -feX Q -fiZ a -pLD a (D b v b ) - <dfin a (58) 
Z (a) = ~QZ a -\X a + D a {D b h b )-\Q 2 h a -\iih a . (59) 

Note that these equations differ from Eqs. (49) and (50) for the comoving Lagrangean u- frame in Ellis 
and Bruni p"7j ] by the presence of the second-order derivative terms as well as the terms algebraic in n a . 
Of importance in the following will also be the first-order evolution equations for the spatial divergences 
of both v a and h a , given by 

(D a v a Y = -lGD a v a -D a h a (60) 
(D a h a Y = {a-l)QD a h a + al^D a v a + aD a {C 1 ) a , (61) 

respectively. 

In terms of X a and Z a the constraint equations ( |l9| ) - ( p2[ ) , (^9|) and ( [43| ) now assume the first-order 
forms 



FLRW -linearised constraint equations: 

= (Ci) a = Z a -\y,v a (62) 

= (C 2 ) a = e abc D b n c (63) 

= (C 3 ) Qb = E ab -D {a n^ (64) 

= (C 4 ) a = D b E ab -\X a + \Q^iv a (65) 

= (C 5 ) a = D a {D b n b )-\X a + \Q i iv a -\Q 2 h a + nh a (66) 

= (C 6 ) a = D a {D b v b ) - |D fc D (a ^ - i9 2 / + p Q . (67) 



The momentum density of the moving matter q a — fiv a induces a spatial gradient Z a in the expansion 
rate of the normals n by (Ci) a . Also note that in the given situation E ab = D/ n 6 \, provided {Cz) ab 
is satisfied and preserved along n; to first order the electric Weyl curvature of (Ai, g, n) is completely 
determined by the distortion of n a . This relation reflects the connection between the tidal field and 
the gravitational scalar potential in the proper Newtonian theory [12]. We remark that in first-order 



approximation the 3-Cotton-York tensor of Eq. (B2) is zero. 

We now find that the set of constraints ( |62] ) - (|67|) evolves to first order according to 

(<7i)<°> = -e(dr + (C 5 ) Q (68) 

(C 2 )^ = (a~l)e(C 2 ) a (69) 

(C 3 ) (afc> = -e(C3) Qb -aD <Q (Ci) h> -(a + i)[eD (a n b > + | A/ D <a z; h > ] (70) 

(c 4 )< a) - -|e(c 4 ) a + i M (c 1 ) a + i/i(c 6 ) a (7i) 

(C 5 )< a > = {a-l)Q{C 5 ) a + aD a {C A )-a\Q 2 {C 1 ) a + (a + \)^{C 1 ) a 

+ f (a + i) /i [ D" (D t t- 6 ) + i 9 /J" 1 X" - | 6 2 / + /i / ] (72) 

(C 6 ) (a> = -e(C 6 ) Q ~|D b (C 3 ) ab + |(C4) a -(C 5 ) a (73) 

(C A ) = -lQ(C A )+D a (C 5 ) a , (74) 



where we have introduced the auxiliary constraint (C A ) ■— D a (Ci) a in order to bring the evolution system 
( |68| ) - (JTJ) into a form suitable for application of the Cauchy-Kowalewskaya theorem on existence and 
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uniqueness of solutions for sets of analytic initial data (see, e.g., Ref. [Q ). As can be immediately seen 



from Eqs. ( |7C| ) and (|72|), we obtain the remarkable result that the set of constraints m 6 \ ) - (61) evolves 
consistently (without restrictions on the background model) {or the special parameter value of a = — g 
in Eq. $3b\). As mentioned before, this is the situation when the dimensionless comoving lapse function 
(SN) is covariantly constant along n. 

If, however, a ^ — |, further constraints arise, namely 

= (C 7 ) Qb := eD^ + f^DV* (75) 



- (C 8 ) ab 
= (C 9 ) a 



i, 



= D (a n 6) +|9D<V> (76) 
= [ (Co) - 3 R } D ( V> , (77) 



the first arising from Eq. (|70|), and the latter two are each the condition resulting from consistency of 
the previous one. Indeed, these propagate along n according to 

(C 7 ) <ab) = (a-l)e{C 7 ) ab + aQD {a (C 1 ) b) -2fi(C s ) ab (78) 
(C 8 )< ab > = ( a -l)9(C 8 ) Qb + aD< a (C 1 ) b > -f (a-i)(C 9 ) Qb (79) 
(C 9 )< Qb > = -|8(C 9 ) ab -i[ (C G )- 3 R](C s ) ab . (80) 

Finally, one can show that for a ^ — i the non- vanishing source term in the square bracket on the 



right-hand side of Eq. (72) can be expressed in terms of the other constraints, provided that 3 R = 0: 

(Ci ) a := D a (D b v b ) + ±e f i- 1 I>y-Z@ 2 v a + tiv a (81) 
= (C 6 ) a + 9 tr 1 [D b (C 3 ) ab - (C* 4 ) Q +D b (C 8 ) ab } - f [ (C G ) - 3 i?]D b D<V> . (82) 



Hence, for u ^ and 0^0 (and assuming that (Co) is satisfied and preserved along n), Eqs. ( |75[) - 
((77]) and (|6J) demand for a ^ — | that 

= 3 i?D (a « b) =>■ = 3 i?£; a6 ; (83) 

with this choice of evolution for N a complete and self-consistent set of variables and equations describing 
linearised scalar perturbations of FLRW dust cosmologies from a non-comoving point of view requires 
either 

• a spatially flat background time slicing, 3 R = 0, or 

• a peculiar velocity field v a with (to first order) vanishing shear, T)( a Vb) = 0. However, this is just 
the conformally flat case which we will further discuss below; the inhomogeneities are pure gauge 
in this case. 

Thus, the conclusion is that there are two consistent inhomogeneous FLRW -linearised cases: 3 R 7^ and 
a = — 5, or 3 R = and no restriction on a. In each case the set of resulting equations is then consistent, 
but in the latter case there are the two extra constraints ( f7^ ) and ( f7^ j that must be satisfied, as well as 
Eqs. MM) - mPK), which must be satisfied in all cases. 



6.2.3 Evolution of spatially inhomogeneous deviations 

It is now convenient to define scaled covariant and gauge-invariant perturbation variables by Jl~7| , ^| 

va 

S a := , Z a := SZ a , (84) 

the first being the dimensionless comoving (with the Eulerian observers) fractional total energy density 
gradient, the second the comoving rate of expansion gradient. The purely scalar content of these variables 
is obtained by taking their magnitudes, or, essentially equivalently, by constructing their comoving spatial 
divergences (l8[ ||: 

5 := SD a S a , Z := SD a Z a . (85) 
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In terms of the latter, Eqs. (p8|) and (^) convert into 

6 = -Z~STJ a [SD a (D b v b )} -eSD a (Sh a ) (86) 
i = -lGZ-±fiS + SD a [SD a {D b h b )]-±e 2 SD a {Sn a )-±[iSD a {Sh a ) , (87) 

which, finally, after some more work, can be combined into a single second-order linear ODE for 6 given 

by 

<5+§6<S-±m<5 = -a|6 j uS'D Q (5w a )-i[2(a-i)e 2 + (C G )- 3 i?] SD a (SYi a ) 

- ae5D a [5(d) a ] . 



In obtaining this relation, use was made of Eq. (30) and the fact that, with Eq. (|37j), we have to first 
order 

[SD a (Sh a ) Y = a[ 6 SD a (Sn a ) + f nSD a (Sv a ) + SD a [S(d) a ] } . (89) 

Equation ( |8S| ) needs to be compared to results given in the standard literature (see, e.g., Refs. || 
and [|6| ), and the Bertschinger relation (|57|) given above. With (Ci)° and (Cq) satisfied initially and 
preserved along n, extra velocity-induced terms occur. The most complex extra source term on the right- 
hand side will arise when a — — | which, as derived above, is the only case that allows for the possibility 
3 R ^ 0. However, even when 3 R = 0, unless a = we have to account in particular for the first source 
term on the right-hand side of Eq. ( pq ) , which directly derives from the momentum density of the moving 
matter. 

As Roy Maartens |3lj kindly points out to us, in the special case when a = (corresponding to 
N = 0) Eq. pH ) can be brought into standard form ( e.g., Eq. (11.1) in Rcf. Q ) by using = fj, — | 6 2 
and defining a new scalar perturbation variable Y := 6 + 2 SD a (Sn a ); the second term then is covariantly 
constant along n by Eq. (^9|). In this sense, computational simplicity singles out the value a = amongst 
all cases with a spatially flat background. The question is as to the physical content of this feature, since, 
unless 3 R 2^0, the value of a can be chosen completely arbitrarily. Solutions to the modified evolution 
equation ( p8[) for 6 can probably only be obtained by numerical means. 

6.2.4 Tilted FLRW models 

We now return to the special case considered above, where D^v^ = 0. This is just a FLRW model seen 
from a non-comoving frame. Indeed, from the non-comoving perspective, an Eulerian observer perceives 
spatial inhomogeneity even if the dust matter flow is exactly FLRW, i.e., if 

= u a = a ab = £u a =► = E ab = H ab , 

because this is then a tilted spatially homogeneous cosmological model . To first order in the deviations 
we find that if this tilted situation is given, then with fi ^ and O ^ it is characterised by: 

(i) from Eq. @ 

= D {a v b) , (90) 

(this is the condition that the spacetime is conformally flat to linear order); and 

(ii) from Eq. (|65|), combined with Eqs. ( |6^ ) and (|30|), 

o = (c 4 ) = -iD a [ 3 i?-(c G )]-|e(dr ; (9i) 

with (C\) a and (Cq) satisfied, the intrinsic 3-Ricci curvature of the time slicing is found to be spatially 
homogeneous. From Eq. ( |67| ) we then have = D a (D b v b ) + \ 3 Rv a . 

6.3 Derivation of the limiting Newtonian-like gravitational field equation 



How does the difference between Eq. fl88|) and the Newtonian form (57) given by Bertschinger arise? The 
key role in standard treatments of FLRW-linearised matter and geometry perturbations is played by the 
PFE, the central equation in the Newtonian theory of gravitation, which, in our context, ties the spatial 
inhomogeneities in the total energy density distribution to the Laplacian of the relativistic acceleration 
scalar potential N. The problem seems to lie in the assumed exactly Newtonian form of this equation 
used by Bertschinger. We suggest that from the point of view of internal consistency this may not be the 
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appropriate Newtonian-like gravitational field equation to use in the cosmological context we consider 
here; that rather this equation should be the PFE modified by inclusion of a momentum density induced 
term. 

Thus, the issue now is to discuss obtaining a consistent set of quasi-Newtonian dynamical equations 
in a cosmological context from the exact relativistic equations. Where can this come from? 



6.3.1 The Raychaudhuri equation 

The Raychaudhuri equation ( |l5| ) becomes the relativistic generalisation of the PFE in the static case, 
when there is no expansion of the timelike reference congruence. Setting O — in Eq. ([l5]) with u 2 «l 
gives the well-known gravi-static equation of gravitational attraction 

D a n a + (h a h a ) = N- 1 D a D a N = \ (i , 

- a direct generalisation of the Newtonian gravitational field equation.^] The problem is that this 
does not work in the generic cosmological context: in a non-static cosmological spacetime this equation 
determines the rate of change of expansion of the model, with the spatial divergence of the relativistic 
acceleration as one of the source terms, rather than directly relating that divergence to the matter energy 
density. 

The Raychaudhuri equation would become the relativistic generalisation of the PFE if we could 
suitably subtract out the rate of expansion, leaving a relation between energy density variation and the 
relativistic acceleration. The problem is that the background expansion is spatially homogeneous but 
the real expansion O is not, so subtracting off the background expansion terms is a gauge-dependent 
procedure. This will only give a gauge- invariant result if O is spatially homogeneous, i.e., if Z a = 0. To 
make this work, we would have to subtract off the expansions terms not of a background model but of 
the real model; but they are spatially varying, and so would add extra terms to the PFE-like equation 
obtained. The resultant extra terms may be negligible in specific contexts (e.g. gravitational lensing 
models), but then that circumstance needs to be shown. 



6.3.2 The Friedmann equation 

The Friedmann equation for our quasi-Newtonian models takes the form 

2 fj, = 3 R + | e 2 . 

If the background FLRW model is spatially flat, then with the line element in the form (|53| ) this relation 
becomes 

2 M = 52 [ + ®\w + *l« ] + I 92 • 

Hence, again if we could subtract off the expansion term | 2 we would have an equation of the form 
we want. However, again that will not work in the generic physically relevant case when O is spatially 
inhomogenous; we can subtract off the background expansion consistently in a position- independent 
way, but not the real expansion, for that is position-dependent. Again, this spatial dependence may 
be negligible in particular circumstances, but that needs to be shown (note that we are carrying out 
a consistent linearised study of the equations; it is not true that linearisation by itself implies we can 
neglect these terms). 

What we really need to do is to take the spatial gradient of this equation to get a covariant and 
gauge-invariant equation of the form we want. However, we already have equations of that kind at hand: 
namely, the D;,i? ab -constraint (Bianchi identity) and the compatibility condition (Cs) a binding the former 
to the shearfree slicing condition (C^)" 6 ; we consider this next. 



6.3.3 The D b £ afc -constraint (Bianchi identity) 

In the case of a spatially flat background, the compatibility condition (65)° between the D;,i? ab -constraint 
(Bianchi identity) (64)" and the zero-shear condition (C3) ab assumes from Eq. (66) the form 

D a {D b h b ) = \X a -\Q^v a . 

25 But generically with active gravitational mass (p, + 3p) (see, e.g., Whittaker or Ehlers jj| ). 
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In terms of the relativistic acceleration scalar potential N, this relation can be re-written as an effective 
Newtonian (PFE-like) gravitational field equation given by 

D a [ D b D b N - ±N n } = - ±NOfiv a ; 

this is vectorial and third-order because it is covariant and gauge-invariant ,p| and it has an extra velocity 
part essentially resulting from the momentum constraint equation (C\) a = 0. We believe that any other 
approach that consistently derives a scalar second-order PFE-like equation in the current context will 
also lead to such an extra term. 



6.3.4 Neglecting the momentum density 

If (as, e.g., in Ref. |27| ) one uses a version of any of these equations without the term proportional to 
q a = /hi", then one has neglected the first-order momentum density due to the non-relativistic peculiar 
motion of the matter. While this term is small, the implication, via (C\) a — 0, is that one has also 
neglected the spatial gradient Z a of the expansion rate of n; for example, the set (|54|) - (|5^) neglects 
this (0a)-EFEj^| which links the matter momentum density — linear in v — to mixed time-space 



derivatives of $ (cf. Ref. |27f| and subsection A. 5 of the appendix) 



To investigate this, we now inquire into the viability of imposing the extra assumption 

q a = /j,v a =0 => Z a = 0; 

thus, we set q a to zero throughout, implying, as \i ^ 0, that v a vanishes. We here find ourselves in a 
slightly delicate situation, as, e.g, Bertschinger Q neglects q a but not v a itself (part of his dynamical 
set is an evolution equation for v). Note that, to first order, we have, e.g., D^g^ = fi'DuVty and 
D a g Q = M D Q i/\ 

From the gauge problem point of view, our extra assumption amounts to choosing a family of spacelike 
3-surfaces T: {t = const} with timclikc normals n which, simultaneously, are shearfree and have spatially 
constant rate of expansion ||. However, for \x ^ and ^ 0, it now follows that if D^aQb) — 0, then, to 
first order, 

D(a«b) = , 

which, for v a ^ 0, we have recognised to characterise tilted FLRW models. In view of this result 
and the need to obtain consistent Newtonian limits from the exact equations of General Relativity in a 
cosmological context, enforcing the PFE by neglecting the first-order momentum density q a of the matter 
in the (Oa)-EFE seems questionable, even though this term is certainly very small. 

It may be that other approaches, e.g., using 'synchronous' gauge conditions, get around this problem 
in a satisfactory way. However, for example Peebles' approach to the CMBR anisotropies, that is largely 
based on this gauge, in fact mixes different gauges — setting up the CMBR anisotropy investigation in one 
gauge, and using another one to establish the PFE relating density inhomogeneities to the gravitational 
scalar potential (see Ref. 0, pp500-6). This seems a problematic procedure, and does not appear to 
solve the issue raised here. 



7 Conclusion 

Put in loose terms, in this article we have asked the question as to what extent we can squeeze a 
generic dust matter source field, which typically generates non-zero 0, a a b, £j a , E a b and H a b, into a 
Newtonian-like reference frame characterised by the conditions (Q) - (Q). The implied invariant physical 
requirement that a quasi-Newtonian cosmology {M, g, n) excludes gravitational radiation provides a 
very strong restriction on the generality of the resulting cosmological model. It may turn out that this 
constraint proves to be sufficiently strong indeed that no exact spatially inhomogeneous solutions (in 
non-comoving description) exist at all. The complexity of the constraint analysis encountered in section 
^ is certainly a reflection of the restrictiveness of the no-gravitational-radiation condition.^] Indeed, it 

26 Again, its purely scalar content is obtained from the spatial divergence of this relation. 

27 The (Oa)-EFE is equivalent to the D f) cr tli '-constraint (Ricci identity) of the 1 + 3 covariant dynamical equations. 

28 Note that this is different from other investigations of cosmological models without gravitational radiation, in particular 
those with H a b = p^ |, because in those cases the zero-magnetic- Weyl-curvature condition was imposed in the matter- 
comoving frame. 
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appears plausible that after eliminating the non-Newtonian spacetime curvature source terms (isotropic 
and anisotropic pressure) by FLRW-linearisation of the quasi-Newtonian equations, the gravitational 
radiation inducing quantities have been switched off, and the system of dynamical equations becomes 
tractable (and possibly solvable) . If it could be shown that non-trivial solutions exist for the linear case 
but not the exact equations, linearisation instability of this dynamical setting would be shown. 

Whether or not this is true, our analysis indicates problems arising in the Newtonian limit procedures 
in the literature based on quasi-Newtonian type conditions of the kind we discuss; use of mixed gauges 
as in Peebles' discussion |37j does not clarify these issues. 

(1) The approximations made in the standard literature, ultimately leading to the PFE of the New- 
tonian theory of gravitation, cannot be regarded as providing an example of a we/Z-dcfincd (intrinsically 
consistent) Newtonian limit of the (cosmological) EFE, as demanded by Ehlers in, e.g., Refs. Q and |J. 

(2) If we use a covariant and gauge-invariant approach, this gives a setting where the ten general 
relativistic gravitational potentials reduce effectively to one, as in Newtonian theory (hence the common 
use of this approach in linearised studies). However, we still have unresolved the issue of what determines 
a suitable choice for time evolution of the N(x l ) — we will get restricted results if we impose a choice 
for N; but we have been unable to see how the integrability conditions generically determine the choice 
of N(x z ) (which is what we expected). Thus, we do not have a solution in sight that satisfactorily 
resolves the problem. This correponds to the unresolved question of how one obtains time-evolution for 
the Newtonian gravitational scalar potential with realistic cosmological boundary conditions. 

(3) Factoring out the subset of (inhomogeneously perceived) spatially flat FLRW solutions in the 
linearised case characterised by the condition T)[ a Vb) = => o a \, = 0,0 we obtain extra velocity-dependent 
terms in the Newtonian-like gravitational field equations resulting from consistent linearisation of the 
relativistic equations, rather than just the PFE that many prefer to use. 

(4) Underlying all this is the difficult issue of coarse graining: in essence we try to describe two 
distance scales (sub-Hubble-radius scales vs cosmological/spatial curvature relevant scales) by just one 
set of (mathematically consistent) dynamical equations. The relation between these different scales of 
analysis (and the corresponding effect on the EFE) needs to be clearly brought out and clarified. 

Astrophysically oriented readers may immediately point out that our use of geometrised physical units 
disguises the true orders of magnitudes of the physical effects induced by various terms in our dynamical 
equations. For example, upon re-establishing the fundamental constants c and 8irG/c 2 (in SI units we 



have 8ttG/c 4 = 2.076 x 10" 43 s 2 /(kgm) ), the constraint (d) a of Eq. (|62|) reads 

= D a 9 - pv a ^ 

where p denotes the mass density of the matter constituents. Similarly, non-standard source terms 
proportional to v a in each of Eqs. (|6^) and ( |88| ) become quite small for non-relativistic magnitudes of v a . 
Thus, in particular cases such as gravitational lensing settings, the extra terms may be negligible. Such 
arguments must, however, be pursued with care; as J L Synge has pointed out (see Ref. ppl76-8), 
if you too easily drop pressure gradient terms on the basis of similar arguments, you may deduce that an 
ocean liner will sink to the bottom of the sea. The argument needs to be carefully made in each context, 
for example, in Zel'dovich-like gravitational collapse scenarios. What we are pointing out is the necessity 
when doing so to look at the consistency of the whole set of equations, rather than looking at just one or 
two equations and ignoring the rest. 
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The spatial rotation e abc D(,d c = => u) a = is of second-order smallness in the first place. 
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A.l Commutation relations 

The following commutation relations for 1 + 3 covariantly defined scalars / and vectors V a have been 
frequently used within our work: 

h b a(Pbf)' = [D +n„](/)-|eD / (92) 

D[„D 6] V c = E c[a V b] + \nv c v [a V b] + h c[a E b]d V d + \ (v d V d ) \i h c[a v b] (93) 

+ i[ (l-« 2 ) M -^0 2 ] h c[a V b] (94) 

D [Q D b] y" = -LE ab V b -\(v b V b )nv a -±[(l-±v 2 )n-±e 2 ] V a . (95) 

Recall that in general the commutation relation for spatial derivative operators acting on / is given by 
Dr a D(,i/ = e abc uj c f (cf. Ref. (l8) ). However, as for quasi-Newtonian cosmologies ui a (n) = 0, within the 
T: {t — const} we have 

D [o D 6] / = 0. (96) 

A. 2 Spatial derivatives of v a 

The irreducible decomposition of the purely spatial covariant derivative of the peculiar velocity v a is 
given by 

B a v b = B {a v b) + | (B c v c ) h ab + e abc [ \ e cde V d v e ] . (97) 

Now defining 

Pab ■= Kb - v a /vv b /v , q ab := v a /vv b /v - \pab , s ab := e abc v c jv , (98) 

individual derivative terms can be expressed by 

B {a v b) = {2v {a /vp i > c b) v^/v + lq ab q cd +p c {a p%}(D c v d ) (99) 

D Q w a = \v a /v 2 D a v 2 +p a b D a v b (100) 

e abc B b v c = v a /v(s bc V b v c ) + 2 S ^ b v^/v(D b v c ) (101) 

D a v 2 = v a /v{v b /vD b v 2 )+p b a V b v 2 (102) 

v b D b v a = ±v a /v(v b /vV b v 2 )+p a b v c D c v b . (103) 



Remember that the spatial rotation of v a as given in Eq. (101) is algebraically constrained through Eq. 
(§)■ 

A. 3 Kinematical and Weyl curvature variables relating to the dust matter 
flow 

The dust matter flow we discuss is moving with average 4-velocity u; we have u a = 0. In deriving the 
following set of relations we made use of Eqs. dl7|), ( p5j ) and (pif): 

&ab(v) = 7 [ D( a i> fc) + | On (a v b) + | 9w 2 (n a n 6 + | /i a6 ) + n (a i> c D| c |U fc) ] 
+ 7 3 [ 5 ("(a + u (a) Db)« 2 + "(a («&) + "&)) (w c w d D c w d ) ] 

" § T 3 [ (1 - h' 2 ) + D ^ C ] [ 2 n (a l 'b) + n a n b + v a v b ] (104) 
Q(v) = 7 [ (1 - \v 2 )Q + D a v a ] (105) 

^ah(w) = 7 [ D [a«fc] + n [a v c T) ]c] v b] - \ 7 2 (n [a + i> [o ) D b] t> 2 + -y 2 n [a v b] (v c v d D c v d ) } . (106) 
As usual, the vorticity vector £j a is given by 

~ a 1 ~abc ~ 



e abc w 6c . (107) 
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In the non-comoving FLRW-linearised limit the kinematical variables of the dust matter flow are given 

by 



&ab = D( Q V b ) 

e = e + D aV a 

Uab = ; 

recall that by Eq. (|25| ) T)[ a v b ] is of second-order smallness. 

The Weyl curvature variables with respect to the dust matter flow u can be obtained from 



(ab) -36 &ab - ^°(a^b)c 



(108) 

(109) 
(110) 



(111) 



7 2 [ (1 + v 2 ) E ab + 2 (n (o - v ) E b ) c v c + {E c dV c v ) (n a n b + v a /v v b / v + p a b) } (112) 



H ab (u) = -t)( a u b) +e cd{a tfa d b) 

= 2 1 2 e cd(a [n b) (E c e v e )+E c b) ]v d . 

A. 4 1 + 3 orthonormal frame dynamical equations 



(113) 
(114) 



For a full consistency analysis, we need the dynamical equations describing quasi-Newtonian cosmologies 
with dust matter source field expressed relative to a 1 + 3 orthonormal frame { e a } ( see, e.g., Refs. |j~tf 
and |Q ). After aligning the timelike frame vector field eo with our reference congruence n, i.e., eo = n, 
one can choose a spatial frame { e Q } according to convenience. Making use of this freedom, the following 
set of quasi-Newtonian relations, which we derive from the general equations given in Ref. p3| by simply 
deleting terms, can be further simplified. 



Commutation relations: 



[e 0) e a J = n a e 
[e a ,ep] = [ 2 ct[ a (Taj + e a ps n &1 



e /3 



Einstein field equations: 

e (6) 






- § 6 2 + (e a + h a - 2 a a ) (h a ) - \ (1 + v 2 ) /* 

(<T <Q e 7 + ri< Q + a< Q ) (h®) + n v {a - *S al3 - e<^ a n \ h s 



= *i? + fe-2 M = (c G ) 

= S a ^e p (Q)~lriv a = (d) Q , 



where 



<Ja/3 

*R 

b a f3 



= e («( a /3>) + b (ai3) - e 7 V (e| 7 | - 2a| 7 |) [np) S ) 
= 2{2e a -3a a ){a a ) - 
= 2 n al n'ji — n 7 7 n a p 



- b a 

2 



Jacobi identities: 



(115) 
(116) 

(117) 
(118) 
(119) 
(120) 



(121) 
(122) 
(123) 



(124) 



e (a a ) = -l(6 aP e,3 + h a +a a ){0) + ±e al3 ~> (ep + h p -2a f3 ){n~ / ) 
eo(n a/3 ) = -16^ + (6^ a e 7 + h^) (fl^) - 6 a/3 (e 7 +n 7 ) (fi'') - 2e*n« Oj (125) 



= e Q ^(e /3 -a /3 )(n 7 )-n Q /3 n' 3 = (C 2) 



= (e f3 -2a )(n a V)+e a ^e (a 7 ) 
Electric Weyl curvature: 

E a p = (e( Q +h( a +a {a ) (hp)) + ^v {a v p) - e 1 \ a n 0)l h s + (C 3 ) Q/3 

= *S a f3 — ^/ID(a^) — (Cg)q/3 ■ 



(126) 
(127) 



(128) 
(129) 
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Contracted Bianchi identities: 



e (fi) = - (1 + ±v 2 )6^ - v a e a (fi) - fi(e a + 2h a - 2a a )(v a ) 
e (v a ) = -\{l-v 2 )Qv a -v l3 {ep-np-ap){v a )-v 2 a a -h' 
+ e a ^ [ n p5 v 5 v 1 + Slpv 1 ] . 



(130) 



(131) 



A. 5 Line element and linearised Einstein field equations for scalar-perturbed 
spatially flat FLRW dust 

Many perturbation calculations use the following quasi-Newtonian line element, where rj denotes a di- 
mensionless conformal time coordinate and x, y and z are Cartesian coordinates comoving with n: 



Infinitesimal line element: 

ds 2 = S 2 [ - (1 + 2$) dr, 2 + (1 - 2*) ( dx 2 + dy 2 + dz 2 ) ] . 



(132) 



The associated commutation functions and components of the Einstein curvature tensor, as well as the 
first-order components of the energy-momentum-stress tensor are: 



Commutation functions: 



" s 

= cr a p = uj a = Q a = n a p. 
Einstein curvature tensor: 



6 = l 



(1 - *) ^ - *|„ 



(133) 



Goo — 

Goa = 
C a a 

G a (3 = 



1 

S 2 
2 

S 2 



3(1-2$) ^ 2 -6|^ + 2*l« |c 



S 



1 77a 



(134) 
(135) 



1 



1 



2 (1-2$) :±22-(i-2*) 



Si, 



no summation, a ^ [3 ^ 7 



2 -f ($ + 2*)|, - 2 - (* - *) m - ($ - *) |77 



52 ( $ -*)|«/3 



Energy-momentum-stress tensor: 

Too = /i 



T 0a = /x w Q 



T Q/3 = 



(136) 
(137) 

(138) 
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